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We analyze a mathematical model of the epizootic of Hantavirus in mice population, including the 
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m : effect of species which compete with the host. We show that the existence of the second species 

■ has an important consequence for the prevalence of the infectious agent in the host. When the 

two mice species survive in the ecosystem, the competitive pressure of the second species may 
lead to reduction or complete elimination of the prevalence of infection. The transition between 
the disappearance of the infection and its presence occurs at a critical value of the competitor's 
population, resembling a second order phase transition in a statistical system. The results provide 
a rigorous framework for the study of the impact of biodiversity in the propagation of infectious 
diseases, and further lends itself to future experimental verification. 
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I. INTRODUCTION: THE ROLE OF BIODIVERSITY 
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t^J- \ In 1993 an outbreak of a severe disease now known as Hantavirus Pulmonary Syndrome (HPS) struck the region 
of the Four Corners, in the North American Southwest, with a mortality in excess of 50% of those affected. Shortly 
afterwards, Sin Nombre Virus (Bunyaviridae: Hantavirus, SNV), the first Hantavirus to be discovered in the Americas, 
• i— i . was identified as the infectious agent responsible for HPS. The main host of the virus was identified as one of the 
most common mammals in North America, the deer mouse, Peromyscus maniculatus. Since then, continuous efforts 
q-i in different areas of science have been made to understand this epizootic, with the ultimate goal of correctly assessing 
the risk to humans. Across all of the Americas, new discoveries of Hantavirus es (half of them being human pathogens 
and some of them responsible for a high mortality ijMills and Childsl 119981) 1 have led to a renewed interest in the 
natural history of the host rodents. 



A simple mathematical model for the spatio-t emporal patterns in the spread of this epizootic has b een proposed 
and analyzed by the second author and others ((Abramson and Kenkrel [2002t lAbramson et all 12003^1 . The model 
takes into account several peculiarities of the Hantavirus-rodent association. An example being, the fact that the 
infection does not produce any known disease in the mice, and consequently does not affect their death rate. The 
role of the sp atio-temporal patt erns of the environment on the prevalence of the infection, which has been found in 
field studies llYates et all 12002ft . has also been taken into account. We refer the reader to Jib ramson and Kenkrel 
200ilAb ramson et all l2003[l for details. That model was able to successfully explain several field observations as 
environmentally controlled phase transitions, thus providing an analytical support to biological hypotheses such as 
the trophic cascade discussed in (jYates et all 12002^) . Among the consequences of the mathematical model that have 
correlations with the field observations, we mention the sporadic disappearance of the infection and the formation of 
"refugia," whence the infection spreads — when conditions change — in the form of waves. 



The model of Abramson and Kenkre (henceforth mentioned as the AK model) predicts a critical transition from a 
state with a positive prevalence of the infection in the population (an infected phase), to a state without infection. 
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FIG. 1 Mean density of P. maniculatus at two sites near Zuni, New Mexico (after ((Yates et all I2002T) . Fig. 7). The rodent 
densities, in animals per hectare, are shown in log scale to emphasize the behavior at low densities. Shown: total population 
(line-dot), infected with SNV (line), critical density (dashed). The time axis counts months from December 1994. See the text 
for a discussion of the dynamics. 



This latter state occurs when a parameter characterizing the environment and controlling the population size drops 
below a certain threshold value, which we call the "critical carrying capacity." The population value associated with 
this critical carrying capacity, m c , is a population threshold for the system: the total population, m, needs to be 
greater than to c , to be able to sustain infection at all. In the terminology of the basic reproduction number, R$, this 
critical condition corresponds to Rq = 1. Above the critical carrying capacity, the infection is able to sustain itself 
(equivalent to the condition R > 1). This threshold phenomenon ha s been documented in the Four Corners region, 
as shown in Fig.^ These population densities have been published bv lYates et alJ l|2002|) . where they discuss the role 
of El Nino on the prevalence of the SNV infection. We plot here the populations on a logarithmic scale to emphasize 
the dynamics for low density values. The line with circles shows the total P. maniculatus density (mean density of 
two nearby sites of the field study), which persists during the 7 years of observation beginning in December 1994. The 
line without symbols shows the density of mice infected with SNV. The dashed line represents the critical population 
able to sustain a positive prevalence of infection, as predicted by the AK model, by us ing approximate parameters 
obtained from the time series (details of the calculation can be found in lAbramsonl <|2004l) ). Let us briefly analyze the 
implications of this graph. At t < 15 the population is above the critical threshold level and, correspondingly, there is 
a positive infected phase. Within this regime, it is conceivable that the environmental carrying capacity has some time 
dependence, but that it remains above the critical value. Around t = 15, the carrying capacity then drops to some 
value below criticality, and consequently the total population m(t) drops as well, approaching an equilibrium which is 
now below the critical value m c . The drop is not monotonic: there are some discrete steps of decreasing population. 
Concomitantly, as m begins to decrease (the indication that the carrying capacity has gone subcritical), the infected 
population to/ begins to disappear sporadically. A few infected mice may be entering by migration, but it is clear 
that the infection is disappearing from the site. After t = 30 there begins a steady population explosion, indicating 
that the carrying capacity has increased, and it is observed that the population grows beyond m c , indicating that the 
system is again above the critical level. The prev alence of infection is expected to recover. This process takes time, 
however, just as in the AK model [see Fig. 2 in l|Abramson and Kenkrel 12002)]. and not before t — 40 do we see a 
positive m/ again. After this, the population remains above critical, and the infection persists. A brief excursion of 
m below 777 c (at t « 70) might be the beginning of a new extinction event, and indeed mi reaches its lowest values 
since those at t = 40. But shortly after this, the time series ends and the analysis can not be carried further. Observe 
that the drop in m takes place in 2001, a y ear that was part icularly dry in the North American Southwest, as can be 
seen in the precipitation data of Fig. 2 in l|Abramsonl l2004|) . 

A threshold in the animal density necessary to sustain a positive seroprevalence has directl y been observed in 
another Hantavirus system, the Puumala Virus in association with the red bank vole in Belgium l|Escutenaire et all 
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2000). It is conceivable, then, that other processes which limit the population size would play a similar role in the 
control of the infection. In real ecosystems mice share the environment with many others species, competing for 
limited resources with some, and being preyed on by others. Indeed, competition and predation stand out as the main 
interspecific relations for the species under consideration. 

The effect of a predator is certainly akin to that of a competitor: both tend to reduce the popu lation under study. 
The role of p r edatio n on the prevalence of a disease has been analyzed bv lOstfeld and Holti ((2004) , and previously by 
iPacker et alJ (|2003^l . As recognized in the former, other phenomena besides predation (such as limited food supply) 
may be the principal regulators of rodent populations. However, under the hypothesis that predation is the main 
factor, they find that the incidence of the disease decreases with an increasing number of predators. An important 
peculiarity of that model is that the number of predators is not affected by the population of prey. As a consequence, 
the number of predators can be used as a control parameter in the analysis of the dynamics. This assumption is 
reasonable for generalist predators, which are sustained by several prey species. A competitor species, on the other 
hand, will necessarily be affected if the shared resources are limited. The mathematical model, in such a situation, 
requires a feedba c k of e ach of t he populations into the dynamics of the other. Various other details of the model of 
lOstfeld and Holtl l)2004j) and of IPacker et al.l l)2003(l do not apply to the problem under consideration in the present 
paper. Among these we might mention: the recovery from infection, the increased mortality of infected animals, and 
the increased predation of them. Certainly, the model could be tailored for specific situations, and for Hantavirus as 
study of the effect of predation will be a necessary complement to the competition analysis that we propose in the 
present contribution. 

Another interesting phenomenon — also relying upon the biodiversity of an ecosystem and affecting the dynamics 
of zoonoses — has received recent attention, as has its relevance to human health . Of particular interest is the ca se of 
Ly me disease, analyzed by mea ns of empirically based computer simulations bv lVan Buskirk and Ostfeldl l)l995tl and 
bvlSchm idt and Ostfeldl (l200lf) . A conceptual model of the phenomenon termed dilution effect, was presented by 
Ostf eld and Keesind l)2000a|) and further developed in ijOstfeld and KeesingLl2000b|) . The dilution effect consists in the 
fact that the prevalence of an infection, in the vector, of an infectious agent (ixodid ticks and a spirochete bacterium, 
respectively, in the case of Lyme disease) is reduced if the population of their hosts contains a diversity of species. 
Ticks (in their juvenile stages, which are the most relevant for the transmission to humans) are non-specific in their 
feeding habits, and therefore if only one species of possible hosts acts as a competent reservoir for the transmission 
of the infectious agent, the existence of other hosts (which are incompetent reservoirs, being unable to transmit the 
bacterium to a tick feeding on them), reduces the probability of infection of susceptible ticks. As a result of this 
"dilution," while the total population of ticks is not affected by the biodiversity of their hosts, the prevalence of the 
infection is reduced. The mechanism by which this reduction takes place in the tick population is, then, different 
from the one we propose in the present contribution. The prevalence of the infection in the Hantavirus host is reduced 
by the presence of competitors at its own trophic level, the species with which the host shares resources. These 
competitors exert a pressure on the host population as a whole, and the infected sub-population may be reduced and 
even led to extinction (even though the total population persists, as we will show below). In other words, with the 
dilution effect there is a reduction in the infection by means of a reduction of the exposure, v.g. of the chance of 
encounter between the infectious agent and a susceptible host. What it is shown in the present work is a different 
phenomenon. Even within the host population, the relative weight of the pathogen is reduced. We stress the ecological 
importance of this feature, rather than its mathematical justification. 

With the purpose of gaining insight into the system, we aim to understand the role of the different processes 
separately, while keeping the model analytically manageable. For this reason we have chosen to make a number of 
simplifications. As mentioned above, we restrict our discussion to competition interactions only, which arise among 
rodent species for a number of reasons. To mention a few, territory, feeding habits, behaviors, etc., all contribute to 
competitive behavior. In general, a particular Hantavirus has a specific host, a single species that acts as a reservoir. 
This host remains infected and maintains the infection in the population by horizontal transmission only. It is known 
that, in some cases, there may be "alternative host species," that are able to host the virus, but for some reason 
are unable to transmit it successfully. Other mammals that come into contact with an infected individual of the 
host species may become infected. However, non- host animals norm ally represent a dead-end for the virus, which is 
eliminated by the immune response of the animal ijPeters et al.lll98fj) . In summary, these species do not play a direct 
role in the spreading of the disease in the wild, but they exercise, nevertheless, an ecological pressure on the host, 
that eventually may affect the dynamics of the zoonosis. 

In the following section we present a generalization of the AK model in the simplest case, incorporating just one 
competitive species that cannot be infected, and studying its effect on the dynamics of the infection. The theoretical 
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insight obtained may suggest ways for controlling the disease agent and, therefore, its incidence. Our results show 
that competition reduces, in a specific way, the prevalence of the infection. Interestingly, this is in agreement with 
(and provides theoretical support to) a hypothesis that has been recently put to experimental test in populations of 
Z. brevicauda, the host of the Calabazo Hantavirus, in field studies in Panama 1 . Similarly, also in Panama, it has 
been proposed that the maintenance of competitive populations may serve to reduce the risk to human populations 
exposed to O. fulvescens infected with the Choclo Hantavirus 2 . The proposal which has been called a "moat" consists 
of an area surrounding human habitation maintaining a diversity of innocuous species, competing with the hosts of 
the Hantavirus. 

As in the AK model, we study the properties of the system as a function of the parameters that characterize the 
(eventually changing) environmental conditions. These parameters are proportional to the carrying capacities of the 
system (one for each species), and therefore provide a natural choice to analyze the dependence of the preva l ence o f 
the infection as controlled by limiting resources, as suggested by the trophic cascade proposed bv lYates et alJ {2002). 
The role of environmental factors that affect rodent density, suc h as plant b i omass and coverage, as determinants 
in the seroprevalence of Hantavirus, has also been suggested bv iBiggs et alJ ((2000) . The potential impact of this 
relationship with the hu man risk of HPS h as been successfully put to test by the use of remote sensing techniques 
and GIS, as reported bv iBoone et al.l (|2000|) . 



II. INTRA- AND INTERSPECIFIC COMPETITION 



We consider two species of mice, one of which is a reservoir of Hantavirus (the host species), and another which is 
not (the alien species). The host is subdivided into susceptible and infected populations. Let us consider that both 
sp ecies interac t by competing for common resources. A usual description of this situation, as discussed for example 
bv lMavl i)l98l|) . is the following model for the population dynamics of the total host, m, and alien, z, populations: 

dm m , . 

— = (b- c )m- — (m + qz), (la) 

-£ = {p-j)z-~(z + em), (lb) 
dt k 

where, for the host species, b is the birth rate, c is the death rate, K is the carrying capacity in the absence of an 
alien population [z = 0), and q is the influence of the alien population; for the alien species, the analogous parameters 
are /?, 7, re and e respectively. We consider that the interaction is not necessarily symmetric, i.e. q =/= e. One of the 
species can take resources from the other in a quantity greater than what it is losing. 

There are four significant equilibria for the system defined by Eqs. lJT]l. the others involve negative populations. 
Two of them consist in one species surviving over the extinction of the other. Another one is the extinct phase where 
all populations are zero. Finally, there is a coexistence state where both species are present, and in this we will focus 
our interest. When only one species persists, the equilibrium population is either: 

m*=K(b~c), (2) 

or: 

z* = «(/?- 7), (3) 
where the asterisk denotes equilibrium states. If the two species coexist, then: 

, _ K(b-c)-q K (P- 1 ) 



1-qe 
n{fi - 7) - eK(b - c) 
1~ ~q~e ' 



(4a) 
(4b) 



1 Suzan, G. (University of New Mexico), personal communication, and also Suzan ct al. (2004|). In their work, only the host was left at 
specific sites, removing the competing species, with the purpose of observing the change in the infection prevalence with respect to 
control sites, where all the competitors coexist. 

2 Kostcr, F. (Lovelace Sandi'a Health System), personal communication. 
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Th e stability of these equilibria is im mediate with the use of linear stability analysis [see, for example, l)Murra\lll993|) 
or ijGuckenheimer and Holmesl fl983)]. The conclusion of such analysis is as follows. If the intensity of the interacting 
competition is not very high, q < 1 and e < 1, then the coexistence state given by Eq. (@J is stable. On the other 
hand, if the competition is strong, q > 1 and e > 1, bistability occurs: the final state depends on the initial conditions. 
Finall y, if q > 1 and e < 1 (or q < 1 and e > 1), only the strong competitor survives. For a more complete survey see 
l)MavMl98lD or llMurravt ITflfll . 

With the same considerations made in the formulation of the AK model, we introduce an internal classification into 
the host population. The internal states consist of the infected subpopulation, m^, and the susceptible subpopulation, 
m s . It is clear that these states are mutually exclusive, and that their sum recovers the total host population, 
to = 7ji s + 77ii . The evolution equations for these are modeled, also as in AK, according to the following field 
observations: infected individuals are generated by pair-wise interactions between susceptible and infected; there is no 
vertical transmission (from infected mother to offsprings); infected individuals die with the same rate as susceptible 
ones (indeed, the virus does not appear to affect any physical or behavioral parameter of the infected individuals) . It 
is worth mentioning that the lack of vertical transmission is, in the AK model, crucial for the character of the phase 
transition controlled by the — environmentally dependent — carrying capacity. If it is relaxed, even infinitesimally, the 
threshold of the carrying capacity becomes zero, and there is a positive prevalence for all values of K . The resulting 
model is now: 

d7n s m s 

= oto — cm s (to + qz) — am s rrii, (5a) 

dt K 

—— = -cm l — (to + qz) + am s m u (5b) 

dt K 

-77 = (P-l)z- -{z + em), (5c) 
dt k 



where a is the contagion rate. The equation for m — m s + mi is the same as in Eq. (jlbjl . The AK model is recovered 
by setting q — 0, and ignoring the equation for z. The existence of the two subpopulations of to makes each of the 
equilibria corresponding to m / split into two, one of them with a positive prevalence (an infected phase) and one 
without infection. The former, for which the three populations are positive, is: 

to* = b/a, (6) 
, _ K{b-c)-qK{p-j) 



1-qe 
k{(3 - 7) - eK(b - c) 
1-qe 



b/a, (7) 

(8) 



Due to the fact that the coexistence equilibrium of Eq. Q can be expressed as: 

m* =K(b-c)-qz*, (9) 

we can cast the fraction of infected mice, \* = m i/ m *> as 

b/a K AK 
X* = 1 213. = i ±£ (10) 

X K(b-c)-qz* K-^z* 1 1 ' 

where K A — b/a(b — c) is the critical carrying capacity of the AK model, below which the infected phase disappears. 
This last result coincides with the AK model when z = 0, so that if K > , and then x* > 0, the infected phase is 
stable, and if K < it is unstable. In the presence of the alien population (z > 0), we can generalize the critical 

parameter as: 

K c = K* K + -S-z. (11) 
b — c 



Let us analyze the relevance of Eqs. i|10[l and Ijlll) . The former, Eq. i|10|) . states that the fraction of infected mice 
in the host population is reduced by the presence of alien mice, the host's competitors. A graphical representation of 
this result is shown in Fig. [21 where we plot x* as a function of the alien population, for several values of the carrying 
capacity K. Fig. [21 contains one of the fundamental results of the present work, clearly showing that the competitive 
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FIG. 2 The fraction of infected mice in equilibrium, x*, as a function of the population of aliens, z (rescaled by q/[b — c]). We 
see that the infected phase is forced to extinction by a critical value of the alien population, that depends on the environmental 
carrying capacity for the host. 



pressure exerted by the alien population has the effect of reducing the prevalence of infection in the host. In fact, from 
Eq 11UI) it can be seen that the prevalence is maximum when no alien population is present, x*i z > 0) < x*{ z = 0), 
and that increasing the amount of aliens will always reduce the prevalence, since dx/dz < 0. A further important 
result of the present calculation is the existence of a critical amount of aliens, a threshold level in the population of 
competitors, that drives the system completely to a non-infected state. This value is given by 

z c = Q/f(6 " c) ' 6 . (12) 
aq 

When the environmental parameter K has a value greater than the critical K c the system has a positive prevalence 
of infection. On the other hand, when K < K c , the competitor population results greater than the minimum 
necessary to force the infected subpopulation to extinction. The point K = K c constitutes a critical point for the 
system, separating two behaviors that qualitatively differ in the stability of the equilibrium of the infected population. 
Observe that time- varying environmental conditions, as well as a heterogeneous landscape, may provide a framework 
where a wealth of spatio-temporal patterns of infection would occur. 

In Fig. |2| we show a phase diagram that provides a complementary picture to that given by Fig. The shaded 
contour plot shows the prevalence (as scales of grey) in the parameter space defined by the two carrying capacities, 
K and n. The region painted white, that includes the origin, is the region free of infection, where the state with 
m* = is stable. The first contour defines the transition line, beyond which a positive and stable prevalence is to 
be found, in a potentially infinite region in parameter space. Two sections of the phase plot are also shown, along 
the lines perpendicular to the axes seen in the contour plot. Each one of the sections represents a projection of the 
prevalence along one of the carrying capacities. Both of them are in the form of a transcritical bifurcation, in which 
two equilibria of different character exchange their stability. The section shown at the top of the figure represents 
the transition as a function of the carrying capacity of the host, K . There is a threshold above which the population 
shows a positive prevalence. A similar picture characterized the AK model (as a function of K only) [see Fig. 1 in 
l|Abramson and KenkreL f2002l) ]. The presence of the competitor reduces the seroprevalence in such a way that the 
bifurcation is not a linear relation between the equilibrium solution and the control parameter. The section shown 
on the right panel represents the effect of the carrying capacity of the competitor, k. In this case, a greater carrying 
capacity inhibits the infection, because it allows a greater population of competitors. 

It is important to remark that the "strength" of the competition with the alien population, q, is the same for both 
susceptible and infected hosts [Eq. (J5J]. Indeed, both subpopulations are reduced as a result of the competition, 
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FIG. 3 Phase diagram of the infection prevalence, in the space defined by the host carrying capacity K and the alien carrying 
capacity k. The contour plot shows the prevalence as shades of grey, with the darkest shades corresponding to higher prevalence. 
The white region is the region of zero infection. The crossing lines indicate sections of the plot, shown in the upper and right side 
plots, where the prevalence is shown as a function of the relevant control parameter. The right panel represents the extinction 
of the infection controlled by the alien carrying capacity. The parameters, in this example, are: b = 1, c = 0.6, q = 0.2, j3 = 1, 
7 = 0.5, e = 0.1, o = 0.1. 

but the infected population suffers the consequences in a stronger manner, becoming extinct in a critical way. This 
phenomenon, analogous to a second order phase transition in a physical system, occurs at a finite value of the density 
of alien mice, z c as shown in Eq. i|12f) . 

The origin of the differential effect by which susceptible and infected hosts feel the competitive pressure differently, 
deserves some discussion. This effect is caused by an important difference between susceptible and infected mice that — 
even though it was already mentioned when the classification was introduced — until now had its true relevance not 
exposed. The difference resides in the ways in which new members of each type are incorporated in the population. 
While the susceptible are born, the infected appear only by contagion of already born (and susceptible) mice. In 
biological terms, there is no vertical transmission, a fact that produces a strong asymmetry between susceptible and 
infected. The lack of a differential death rate between susceptible and infected animals is also a peculiarity of this 
system, which is uncommon in other epizootics. These two ingredients, a difference and a similitude between the two 
classes of host animals, are responsible for the differential effect. 



III. CONCLUSIONS 



We have studied a model of the epizootic of Hantavirus in a system composed of two rodent populations: a host 
of the virus, which remains chronically infected by horizontal transmission, and a competitor species that does not 
host the virus. This situation is paradigmatic of many Hantavirus-rodent systems, in which a single species hosts the 
vir us in a competent way, and se veral o ther species compete w ith the host. The model generalizes a previous model 
bv lAbramson and Kenkra l)2002f) [also l)Abramson et all 12003)]. that described the dynamics of the host only. The 
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interaction between the two rodent species, since it is a competition for limited resources, produces a mutual limitation 
of the populations. In particular, the second species exerts a regulatory effect on the host. This competitive pressure 
is the same for both susceptible and infected hosts (as it should, since the infection with Hantavirus does not affect any 
biological parameters in the hosts). Nevertheless, a reduction of the seroprevalencc in the host species is observed, 
up to a point where it completely disappears from the system. This phenomenon can be described as a critical 
transition between an infected phase and a phase free of infection, controlled by the — environmentally determined — 
carrying capacities of the two species. In the space defined by these carrying capacities spatially or temporally 
varying environmental conditions would be described by a path, which might cross the transition line, thus dr i ving a 
corr esponding change in inf ection phase. This is analogou s to the resu l ts disc ussed bv lAbramson and Kenkrd |2002) 
and I Abramson et al.1 l)2003l) . as well as, conceptually, bv I Yates et al.l l|2002() . Furthermore it is interesting that a 
change in the carrying capacity of the competitor only (a vertical path in the phase space of Fig. |3J) is enough to drive 
the infection to extinction. In a real world situation, where the carrying capacities exhibit temporal as well as spatial 
heterogeneity, a wealth of dynamical transitions are to be expected. 

Alternatively, the competitor's population can be used to define an effective carrying capacity of the host, which 
exhibits the same threshold behavior. Our model has been intentionally kept simple to allow a direct analysis of the 
consequences by analytical means. This thus provides a manageable framework to study a necessary step towards the 
modeling of realistic systems, which certainly contain several species of rodents and small mammals in competitive 
interaction, as well as other species in the t rophic web. T he study of similar systems, for example: when competing 
species are organized in a hierarchical way iTilmantllQQl : the role of fluctuations when populations are small; etc., 
are currently under study and will be presented in forthcoming contributions. 
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